DUO TO SMOOTHING EFFECT, NEGLECTING BOTH OF LOCALLY HIGH VALUES AND VARIANCE OF Geochemical DATA, VISUALIZATION OF Geochemical ANOMALIES VIA COMMON MOVING AVERAGE INTERPOLATION APPROACHES (E.G., IDW, KRIGING AND ETC.) IS PROBLEMATIC. IN ADDITION, BECAUSE OF MULTIFRACTAL NATURE OF MINERALIZATION, THE QUANTITY OF Geochemical ENRICHMENT AND/OR DEPLETION, CAN BE EXPRESSED VIA A MULTIFRACTAL VARIABLE, CALLED SINGULARITY EXPONENT. ONE MAY SYNTHESIZE THE SINGULARITY EXPONENT WITH MOVING AVERAGE INTERPOLATION APPROACHES TO ADDRESS THE DRAWBACKS OF COMMON CONTOUR MAPPING METHODS WITH MULTIFRACTAL MOVING AVERAGE INTERPOLATION. HOWEVER, IMPLEMENTATION OF MULTIFRACTAL MOVING AVERAGE INTERPOLATION METHODS IS CHALLENGING. IN THIS PAPER, AN EASY-TO-USE ALGORITHM FOR EXECUTION OF MULTIFRACTAL MOVING AVERAGE INTERPOLATION METHODS WAS PROPOSED. FOR ILLUSTRATING THE ADVANTAGES OF THE PROPOSED ALGORITHM, WE USED A STREAM SEDIMENT DATA SET OF AHAR DISTRICT. RESULTS REVEALED THE EFFICIENCY AND SUPERIORITY OF THE APPLIED MULTIFRACTAL INTERPOLATION USING THE PROPOSED METHOD OVER THE ORDINARY APPROACHES.